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Abstract 

We study the dynamics of a Fermi gas with a Coulomb interaction potential, and show 
that, in a mean-field limiting regime, the dynamics is described by the Hartree-Fock equa- 
tion. This extends previous work of Bardos et al. [3] to the case of unbounded interaction 
potentials. We also express the mean-field limit as a "superhamiltonian" system, and 
state our main result in terms of a Heisenberg-picture dynamics of observables. This is a 
Egorov-type theorem. 

1 Introduction 

The Hartree-Fock equation is a fundamental tool, used throughout physics and chemistry, for 
describing a system consisting of a large number of fermions. Despite its importance for both 
conceptual and numerical applications, many questions surrounding it remain unsolved. One 
area in which significant progress has been made is the microscopic justification of the static 
Hartree-Fock equation, which is known to yield the correct asymptotic ground state energy of 
large atoms and molecules; see [1,7,8,10,12,13]. The time-dependent Hartree-Fock equation, 
which is supposed to describe the dynamics of a large Fermi system, has received less attention. 
To our knowledge, the only work in which this equation is derived from microscopic Hamiltonian 
dynamics is [3]. The Cauchy problem for the time-dependent Hartree-Fock equation has also 
been studied in the literature; see [2,5] and especially [19], where the Cauchy problem is solved 
for singular interaction potentials. 

A key assumption in [3] is that the interaction potential be bounded. A goal of this article 
is to extend the result of [3] to a class of singular interaction potentials, which includes the 
physically relevant Coulomb potential. We also describe how this mean-field result can be 
formulated as a Egorov-type theorem. 

A system of N fermions is described by a wave function ipN-{xi, . . . , Xn) S /\^ L^(R^, dx) 
which is totally antisymmetric in its arguments. The dynamics of tpN is governed by the usual 
Schrodinger equation. In order to obtain a mean-field limit, the Schrodinger equation is rescaled 
with N. In this article we adopt the scaling of [3]. The Schrodinger equation reads 

idt^Nit) = HN^N{t), (1) 
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where the iV-particle Hamiltonian Hj^ is defined by 



N 1 

Hn ^h, + —^w{xi-Xj). (2) 



i<3 



Here, hi is a one-particle Hamiltonian acting on the z'th particle, typically of the form hi = 
— Ai + where A is the three-dimensional Laplacian and v is some external potential; w is 

the interaction potential. Under the assumptions on v and w we make below, it is easy to see 
that i/jv is a well-defined self-adjoint operator with domain /\^ H^(R'^). 

We briefly sketch our main result. Consider a sequence of A'^ orthonormal orbitals (pi, . . . , ip^, 
where ipi is a one-particle wave function. This defines an iV-particlc fermionic state through 
the Slater determinant 

tpN ■= ipi A ■ ■ ■ A(pN ■ 

Let ipN{t) be the solution of the Schrodinger equation ([T|) with initial state ipN- In general, 
tl'Nit) is no longer a Slater determinant for t 0. However, one expects that this holds 
asymptotically for large N: 

■ipN{t) w ipi{t) A • • • A (pN{t) . 
Here the orbitals ipi{t), . . . , (ppf{t) are supposed to solve the Hartree-Fock equation 

N N 



idt(pi 
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Our main result (Theorem 15.31 below) is a precise formulation of this asymptotic behaviour. 

This result is of some physical relevance for studying the dynamics of excited states of 
electrons in large atoms or molecules in the Born-Oppenheimer approximation. Consider a 
molecule consisting of K nuclei at fixed positions . . . , Rk G M'^, as well as N electrons. The 
Hamiltonian is given by 

^ / K 2 iKT \ 2 



^-^\ ' ^-^ \xi — Rk\ j ^-^ \xi 



Here, ca? is the elementary electric charge which we rescale with iV. The electric charge of 
nucleus k is e^Nzk^ where zi, . . . ^zk are constants chosen so that X^feLi ^fc = 1- This means 
that the molecule is electrically neutral. If we choose — eo/v^, for some fixed eo, the 
Hamiltonian becomes 



N / K 



E 



^ \xi — Rk\ N ^ \xi — Xi\ 

The scaling of the elementary electric charge cn may be justified by the fact that the fine 
structure constant a = ef^ = e^/N is small, i.e. TV = 0{a~^). In fact, a « 1/137. 

One problem in the above model, as well as in the works [1,7,8,10,12,13], is that, as N 
becomes large, relativistic effects should be taken into account. Indeed, a simple argument 
shows that the average speed of the innermost electron of an atom with atomic number Z 
behaves like Za (in units where the speed of light c = 1). Another problem in applying the 
time-dependent Hartree-Fock theory to the dynamics of excited states is that the interaction 
with the radiation field is neglected. This interaction is responsible for the relaxation of excited 
states to the ground state of the molecule. 

A physical scenario that is quite different from the large atom or molecule described above 
is an interacting Fermi-gas confined to a box of fixed size. As discussed in [6, 14], the natural 
scaling in this situation may be viewed as a combination of mean-field and semiclassical scalings. 
This problem was first studied in [14, 18]. The authors show that the limiting dynamics is 
governed by the Vlasov equation. These results were somewhat sharpened in [6], where the 
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authors compare the Hamiltonian dynamics with the dynamics of the Hartree equation, and 
derive estimates on the rate of convergence. 

Finally, we outline the key ideas of our proof. It relies on the diagrammatic Schwinger-Dyson 
expansion and Kato smoothing estimates developed in [9]. The main steps are: 

(a) Use the Schwinger-Dyson expansion to express the Hamiltonian time evolution of a p- 
particle observable. 

(b) Show that, in the limit N ^ oo, only the tree terms of the Schwinger-Dyson expansion 
survive. 

(c) Show that the time evolution of a p-particle observable under the Hartree- Fock equation 
converges to the tree terms of the Schwinger-Dyson series as oo. 

Steps (a) and (b) have been addressed in [9]. Thus, the proof below consists in doing step (c). 

The article is organised as follows. In Section [2] we introduce the Hartree-Fock equation, 
discuss its Hamiltonian structure and prove a Schwinger-Dyson series for its time evolution. In 
Section [3] we rewrite the Hartree-Fock equation using density matrices. Section H] is devoted to 
a discussion of the key properties of Slater determinants. After these preparations, we state and 
prove our main result in Section [5l The final Section \E\ is devoted to a Egorov-type theorem, 
which describes the microscopic dynamics as a quantisation of a classical "superhamiltonian" 
theory. 



Conventions 

In the following, the expression "A{t) holds for small times" is understood to mean that there 
is a constant T such that A{t) is true for all \t\ < T. The precise value of T can always be 
inferred from the context. To simplify notation, we assume in the following that t ^ 0. 
The norm of a Hilbert space H is denoted by ||-||. We denote by 

the symmetric/antisymmetric subspaces of the tensor product space 7i®". Here, P± is the or- 
thogonal projector onto the symmetric/antisymmetric subspace. The Banach space of bounded 
operators on H with operator norm is denoted by {B{Ti.), ||-||), and the Banach space of trace- 
class operators on H with trace norm is denoted by {£^{H), 

We use the notation a|^^ e i3(H'^") to denote a p-particle operator a^P^ G B{TL®p) acting 
on the particles ii, . . . ,ip € {1, . . . ,n} in n-particle space. Similarly Tr^j . ^^ denotes a partial 
trace over the degrees of freedom of particles ii, . . . ,ip. 

A time subscript of the form (•)( is always understood to mean time evolution up to time 
t of (•) with respect to the appropriate free dynamics. We shall explain this in greater detail 
whenever this notation is used. 

The symbol C is reserved for a constant whose dependence on some parameters may be 
indicated. The value of C need not be the same from one equation to the next. 



2 The Hartree-Fock equation 

For simplicity of notation, we only consider spinless fcrmions in the following; the one-particle 
Hilbert space is 7i := L^(IR.^,da:) = L^(M'^). Merely cosmetic modifications extend our results 
to the case of spin-s fermions for which the one-particle Hilbert space is (R'^) (g) C^^+^ . To fix 
ideas, we consider the free Hamiltonian h —A and a Coulomb two-body interaction potential 
w{x) := By a simple extension of the results of [9], Section 8, our results remain valid 

for a free Hamiltonian of the form h = —A + v and a two-body interaction potential w, where 
w is even and v,w G L°°(R^) -I- L^(K^) are both real. Here, L^, denotes the weak L'^-space (see 
e.g. [16]). In particular, we may treat Hamiltonians of the form 
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2.1 Some notation 

It is convenient to state the time-dependent Hartree-Fock equation in terms of an infinite 
sequence of orbitals \1/ = (V'i)ieN which is an element of the Hilbert space 

fi := /2(N;L2(R3)) = ?2(N)(g)L2(M3). 
To simphfy notation, we set a — {x,i) and write ^'(a) = tpi{x). Furthermore, we abbreviate 

/ da := / da;, S{a — a') :— 6ii'6{x — x'). 
The scalar product on TL is then given by 



Let a(p) G B{n®P) and define h^P'^ £ B{ii®P) through 

a^^') := l(i2(N))®p ®a'^P^ . 

We have the identity 

||a(P)|| = ||a(P)||. (5) 

Furthermore, 

ii,...,ipeN 

2.2 Hamiltonian formulation of the Hartree-Fock equation 

The time-dependent Hartree-Fock equation for 4* reads 

It is of interest to note that ([7]) is the Hamiltonian equation of motion of a classical Hamiltonian 
system with phase space F := /^(N) ® Y{^{M?). 

Define the map A from closed operators A^^^ on Ti^-* to "polynomial" functions on phase 
space, through 

= / dai • • • dap d^i • • • d/3p ^{up) ■ ■ ■ *(ai)^(P)(ai, . . . , a^; /3i, . . . , /3p) *(/3i) • • • «'(/3p) , 



where v4(p)(ai, . . . , ap; /3i, . . . , /?p) is the distribution kernel of A^p^ (see [9] for details). We 
denote by 21 the hnear hull of functions of the form A(A(p)), with A^p^ € Bifi^^'). 
The Hamilton function is given by 

H := A{h) + iA(W) , (8) 

where 

with 0:2) := ^{x2, xi) and W is the two-particle operator defined by multiplication by 

w{xi — X2). Written out in terms of components, ([8]) reads 
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Using Sobolev-type inequalities, one readily sees that H is well-defined on F. 
A short calculation shows that the Hartree-Fock equation is equivalent to 

The symplectic form on F is given by 

LU = i J da d§(a) A d^'(Q;) , 

which induces the Poisson bracket 

{*(«), 'I'(/3)} = i5(a-/3), vI/(/3)} = *(/3)} - 0. (9) 

Thus, for two observables A,B e 21, 

The Hamiltonian equation of motion on F is the Hartree-Fock equation ([7]) • 

The conservation laws of the Hartree-Fock flow can be understood in terms of symme- 
tries of the Hamiltonian ([5]). One immediately sees that (O is invariant under the rotation 
^ I— > ([/ ® 1l2(r3))5', where U e B{P{N)) is unitary. A one-parameter group of such unitary 
transformations is generated by linear combinations of the functions Re(V'i , V'j) lm{ipi , ipj), 
which Poisson-commute with the Hamiltonian ([5]). By Noether's principle, it follows that 
{ipi , ipj) is (at least formally) conserved. The energy H is of course formally conserved as well. 

In order to solve the Hartree-Fock equation ^ with initial state ^, we rewrite it as an 
integral equation 

V-.W = e-"''^^ - i r ds ^((w; * - {w * . (10) 

The Cauchy-problem for (fTO|) was solved in [19]. We quote the relevant results: 

Lemma 2.1. Let ^' G Til. Then (fTO|) has a unique global solution ^'(•) G C(R; H). Furthermore, 
the quantities {ipi,ipj) are conserved. In particular, ||^(i)|| — \\'^\\- 



2.3 A Schwinger-Dyson expansion for the Hartree-Fock equation 

Our main tool is the Schwinger-Dyson expansion for the flow of the Hartree-Fock equation. 
We use the notation (•)( to denote free time evolution generated by the free Hamiltonian A{h). 
Explicitly, 

^t(Vi,V2,...) = A(e-'*Vi,e-"V2,...)- 

Lemma 2.2. Let A G 21, > 0, and ^'(i) he the solution of pn]) with initial data . Then, 
for small times t, 

A{^{t)) = M^)^ l"ds^{A{W),A^,}{^{s)) 



oo 

fo 2' J A" it) ^ J 



uniformly for "if £ B^, -.^ {"ii E H : i^} . 

Proof. The proof of Lemma 7.1 in [9] applies with virtually no modifications. One uses ([5]), the 
identity 

A(W)t = A{Wt) = A{{Wt{l-E)y), 
and \\E\\ =1. □ 
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3 The density matrix Hartree-Fock equation 

From now on, we only work with orthogonal sequence of orbitals 

/C := {vl/eTi : =0, z^j}. 

By Lemma ^ G IC implies that ^'(t) e K, for all t. To each sequence of orbitals ^' we assign 
a one-particle density matrix 

It is easy to see that this defines a mapping from JC onto the set of density matrices 

V := {jeC^iH) : 7^0} . 

Furthermore, 

ii7*iii = m'- 

Conversely, one may recover ^E* from 7*, up to ordering of the orbitals, by spectral decomposi- 
tion. Furthermore, ^ implies that 

A(5(^'))(*) - Tr(a(P)7|P). (11) 

Let be a solution of ([7|) with initial data 5* and write 

lit) = 7*(t) ■ 

Then a short calculation shows that 

ida = [/i,7]+Tr2[>V,7®7] , (12) 

which is the Hartree-Fock equation for density matrices. As an integral equation in the inter- 
action picture, this reads 

7(t) = e-'*Se"''-i / ds e^'^*-'*)'' Tra [W, 7(5) ® 7(5)] e'(*-^)^ (13) 
"'0 

Sometimes it is convenient to rewrite this using the shorthand 

^(t) := e'*''7(^)e-"^ (14) 

Then (fT^ is equivalent to 

lit) = l-i [ ds Tr2[W„7(s)®7(s)] • (15) 



Lemma 3.1. Let ^(t) be the solution of (fTO| . Then 7,f(f) solves (fT3l 
Proof. Let a^^) = a e BiH). Fr om Lemma we get 



Ai&)i^it)) = Ai~at)i^)+ / ds{A(W),A(ai_,)}(v|/(s)). (16) 

Jo 

Now dH) and imply 

{A(W),A(a)}(*) = iA([W, a® !])(*) 

= iTr([>V,a(g) l]7*(8)7v[,) 

= -iTr((a(g) l)[yV,7M, 7*]) . 

Thus HH) reads 



Tr(a7v[,(t)) = Tr(at7^,)-i / ds Tr((at_, ® 1) [W, 7>j,(,) ® 7*(,)] ) 

Jo 

= Tr(ae'*S*e-'*'^) -i^ ds Tr(a e-'^*"^)'' Tr2 [W, 7*(.) 7*(.)] 



Since a € BiH) was arbitrary, this is equivalent to (fT5|) . □ 
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4 Slater determinants 

The Hartrcc-Fock equation naturally describes the time evolution of quasi- free states [4]. Let 
ujj be the quasi-free state corresponding to the one-particle state 7 G I?. Define 

^'^\xi,...,Xp;yi,...,yp) := u)^ {a* {vp) ■ ■ ■ a* {yi)a{xi) ■ ■ ■ a{xp)) , 

where a*{x), a{x) are the usual creation and annihilation operators of the CAR algebra over H 
(see e.g. [4]). The quasifreeness of w-y means that 

-1^^\xi,...,Xp\yi,...,yp) = A.et{'y{xi\yj))^ . . 

In other words, the operator kernel of 

^(P) = ^^Y^ip) ^ (17) 

where 

For the following calculations it is convenient to introduce the symbol which is equal 

to sgno- if ii,. . . ,ip are disjoint and there is a permutation a G Sp such that {ii,. . . ,ip) = 
(i<7(i)) • • • )i<T(p))) and equal to otherwise. Also, for the remainder of this section, summation 
over any index appearing twice in an equation is implied. 

Lemma 4.1. LetjeV with Tt'j = 1. Then Tij'-p^ < 1. 

Proof. There is an orthonormal basis (<Pi)ieN and a sequence of nonnegative numbers (A,)igN 
such that J2i ^ ^ and ■y = J2i Therefore, 



This yields 



7^^^ = <.'.'.t^n • • • vOi^h • • • VjJ 

Tr7(P) = £■;;::■; • • • X^Ji^k^ ■ ■ ■ Si^kJnk^ ■ ■ ■ Sj^k^ 

disjoint 

^ Ail • • • = 1 • 



□ 



Next, we introduce a special class of quasi-free states, described by Slater determinants. Let 
AT G N and take an orthonormal sequence of orbitals '^n = {<fi, ■ ■ ■ ,<fN)- We define the Slater 
determinant 

Note that the normalization is chosen so that ||S'($jv)|| = 1. The corresponding A^- particle 
density matrix is 

r^V := |5($jv))(5($;v)|. 
One finds for the p-particle marginals 

rj^^ := Trp+i...NrN 

= Trp+i...jv ^eji...»„£ji...j„ \<fn fiN){fji <fjjv\ 

(TV — »)' 

= \vn^---^vip)(vh^---^ I • (18) 
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In order to relate the sequence $Ar to the results of the previous sections, we define the 
normalised sequence 

■■= -j={(pi,...,ipN,0,...). (19) 

V Jy 



Thus, 4* AT G /C and ||'I'Ar|| — 1. It is trivial to check that ^'jv(0 is a solution of ([7|) if and 
only if $Ar(i) is a solution of ([3]). Similarly, '^wit) is a solution of pll)l if and only if <i>Af(t) = 
(931 (i), . . . ,LpN{t)) is a solution of 



i /-^ ^ 

"'O j=l 

Next, from we find 

Thus ini) implies that 

7^^^ - j^sttkn^■■■^^^.){^n^■■■^^^^ - (21) 

Thus, Slater determinants determine quasi-free states by their reduced p-particle marginals. 
The normalisation (^) differs slightly from the usual normalisation 1 of quasi-free states, 
but in the limit N ^ 00 this difference vanishes. Note also that 

117^11 = ^. (22) 

This is a special case of the well-known statement (see e.g. [11]) that || Tr2...Arr|| ^ N^^, for 
any A^-particle density matrix F 0. 

5 Proof of convergence and the mean-field limit 

We now turn to the proof of our main result. We use the graph expansion scheme for the 
Schwinger-Dyson expansion developed in [9]. 

5.1 The Schwinger-Dyson graph expansion 

For the convenience of the reader we summarise the relevant results of the graph expansion 
in [9]. For details and proofs we refer to [9]. Let a'^' G BiTi.'f^) and define its extension 
AAr(a(P-') to iV-particle space Ti^-* through 















1 is 





N <p. 



(23) 



(5($^),A^(a(P))^(a>^)) = Tr(a(f)7^^)) . (24) 
The Schwinger-Dyson series for the time evolution of A7v(a^^^) is given by 

00 A; ^ 

e»^«A^(a(rt)e-'*«« ^ ^ ^ 2. A^(Gf '^(^(p))) . (25) 

k=Q i=a 



^This can also be inferred from H22| l by writing F as a linear combination of projectors. 
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The Hamiltonian i/jv with mean- field scahng is defined in ([2]). The series (|25|) converges in 
norm for all times t; for small times t the convergence is uniform in N. The {p + k — /)-particle 
operator Gj'"'''' (a'-^-') is a shorthand notation for 

Gf''^(a(f)) / diGf/^(a(P)), (26) 

JA''{t) 

where t = (ti, . . . ,tk) and A''(i) is the fc-simplex {{ti, . . . ,tk) : < tk < ■ ■ ■ < h < t}. The 
operators G['''^\a^'^^) are recursively defined by 



'p-\-K — L — l 



iP_ J2 [w^jM^GiTu-ttlli'''^'^^^^^^ (27) 

l^i<j^p+k-l 

as well as G4^°'°^(a(P^) := aj^\ If Z < or ^ > fc then G^^^;;'^ (^(a^^')) = 0. Here, as before, a time 
subscript denotes free time evolution: 

The operator Gp''^(a(p)) may be written as a sum over graph^: 

= I E diG(:^,f5ja(^)), (28) 



Qe3(p.k,l) 



where iq £ {0, 1}, and Ag(<) C [0, i] , and =S(p, /c, ^) is a set of graphs satisfying 

|i?(p,fc,OI < 2''/^') P^ + ^^Vp + fc-Z)'. (29) 



J, 

The operator G^J^''^.^.^^^ (a^P)) is an elementary term, indexed by the graph Q, of the form 
where r = 0, . . . , k. 

The operator norm of g|'^'''' (a'^-*) may now be estimated using the dispersive estimate 

dt |||x|-ie"'^V'f =^ TrllV-lp. (31) 
Going to centre of mass coordinates and using Cauchy-Schwarz, one sees that pip implies 



ds \\W,,^s4 ^ ]f^\M- (32) 

Together with the estimate ()29p . it is now easy to argue, as in [9], that ([25]) converges uniformly 
in A'' for small times t. Moreover, the large- iV asymptotics of the Schwinger-Dyson series 
is given by the tree terms: for small times t we have 

C5C 

e''"-Aj,{a^P^)e~''"- = ^ Ajv(Gl'^°)(a(''))) ijv(i) , (33) 

fe=0 

where ijv (0 1 corresponding to the sum of all "loop terms" , satisfies the estimate 

WL^m ^ G(p,«,t)||a(f)||iV-i, (34) 

for small times t. 



called graph structures in 
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5.2 Convergence of the Hartree-Fock time evolution to the tree terms 

We now give the main argument of our proof. We show that the Hartree-Fock time evolution is 
asymptotically {N oo) given by the tree terms (i.e. the terms ^ = 0) of the Schwinger-Dyson 
series This result is summarised in Lemma [5.21 below. 

Let ^ = {'>Pi)iZi G and denote by ^(t) the solution of the Hartree-Fock equation (fTU]) 
with initial data ^. Let 7(t) = "f^{t) be the associated one-particle density matrix. 

By choosing A — A{d^P^), a^P^ G B{Ti!"^^), in Lemma and mimicking the proof of Lemma 
13.11 one finds that 

Tr(a(rf7W^^) = Tr(ai^)7®'') - i T ds Trfai^J, Tr,+i [W.^+i, 7(5)®'''+'^] ) (35) 

■^0 ^=l ^ ' 

It is convenient to use the representation 7(t) defined in pljl . Using the substitution 
in ((35l) we get 

Tr(a(P) 7(t)®P) = Ti {a'^P^ j^p) - i ds V Trfa^f) Trp+i [W,p+i,„ 7(5)^(^+1'] 

Recall that W^j{l - Eij). Also, Ei^ commutes with W,j and with 7(s)®(p+i). Thus, 

5](p)^(p) implies 



7(i)®f = 7«5P sL^) - i / ds Trp+i 

Jo 

On the other hand, using 



(36) 



.7=1 ^ 



we find 



Trp+i 



^l^.p+i,„7(s)®(P+i)S 



(p+i) 



= Tr 



= Tr 



P+i 



P+i 



^ iy.p+i,„ 7(s)®(P+i) 1-5] S 



1=1 
p 



p 



i=l 



P 



.(P) 



Together with ([55]) this yields 

7(t)®fsL^) = ^®PE(_P)_i / ds Trp+i 
with an error term 



^W.p+i,.,7(s)^(^+^^S 



+1) v(P+i) 



i?p(i), (37) 



Rpit) --^ -i / ds Trp+i [W,p+i,„7(s)®(''+i)£;jp+i 



.(p) 



(38) 



The partial trace is most conveniently computed using operator kernels. We find 
(Wip+i.s7(s)®(P+^)£;jp+i) [xi,..., Xp+i;yi, . . . , yp+i) 

Wsix^, Xp+i;zi, Z2) j{s){zi;y^) ^{s){xj;yp+i) 7(s)(z2; : 



dzi dz2 



n 7(s)(a;r;2/r) 
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so that 



dzi dz2 dz^ 



Wsixt, Z3; zi, Z2) j{s){zi;y,) j{s){xj; Z3) j{s){z2;yj) 



n 7(s)(a;,. ;?/,.) 

= (7j(s)W^y,s7('S)'^^) {xi,...,Xp]yi,...,yp). 
The second term of the commutator in ([55]) is the adjoint of the first and we get 



1 ^ A-LA^ ^ ^ 



(39) 



We proceed to show that, up to an error term, the expansion of the Hartree-Fock time- 
evolution is equal to the tree terms of the microscopic quantum-mechanical evolution. Let 
a(P) e B{l-6^^). Using ^ we find 



Tr 



®p Y-iv) 



Tr (ajP) 7(t)®P eL^^ 

Tr (ajP) 7®'' eL^)) +i [dsY, Tr ([w-,,+i,„ 4''^ ® 1 
Tr (a^^^ Rp{t) 



t P 



i=l 



7(s) 



Iterating this times yields our main expansion 

K-l 



Tr 



{4'^ lit) 



dt Tr 

fc=0 l^i#jXp+fc 



(40) 



where 



4W ^= -i / dtTr(G['iZ..tM'^)Utk+i)W,,,,^,^^^^^^^ 

JA'« + i(t) ^ 

- W,,,,,, 7,(^/0+1) G(:;f:)..,Ja(^))7(t.+i)^(^+^) S^^'^^ 

We now derive a bound on R^j{t). Let us concentrate on the first term, which we rewrite 
using the renaming t^+i — > s as 



dt I ds Tr 

A''(t) Jo 



(Gi:^/"i.,t.(«^^')7,(s)t^.,,.7(s)®<^+'^ 



.(p+fc) 



(41) 



where Ai := min{ii, . . . , tk\. The idea is to use a tree expansion on 7(5). 
Lemma 5.1. Let a^P^ e B(Ti.^f^). For small times we have the tree expansion 

Tr(a('')7(i)®PsL^)) - E / di Tr(T,('=^(a(''))7®(''+'=)EL^)) , (42) 
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where T^^"' is the linear operator defined by T^'^\a^P^) '■— a^^^ and 

p+k-l 



i=l 

Proof. Lemma [2 . 2 1 applied to ^ = A(a'^P') yields 



k=0 



The claim then follows by noting that S^-'o'^^ = pla^^'^ and that J2i=i ^ ^ip+i.t^ commutes 



with The convergence of the series is shown below. 



□ 



Thus (01]) is equal to 



E 



k'=0 



dt 



ds 



dt 



Next, we recall from ([SS)) that G|''t'°'' (^(a^^'-') can be written as a sum over tree graphs 
Q G k, 0) of elementary terms of the form ([50]) . Also, since the definition of t/^';^..,* Jo^p^) 
is the same as the definition of Go!^t°|'..._t^ {'^^^^) with W replaced by W, we immediately get that 
'^tu...,tki^^^'') equal to a sum over tree graphs Q G fc,0) of elementary terms of the form 



This implies that the series (|42|) converges for small times. 

Applying the tree expansion to both c'i'^t^l -t^{a^P'>) and 7(s)®(p+'=) in gl]), we see that (jH 
is equal to 



fc'=0 



E 



Qg^(p,/c,o) Q'e^(p+fc,fc',o) 



dt 



ds 



dt' 



A|',(s) 



Tr 



where A, B, C are operators that depend on (Q, Q', k, k',t,t). A, B and C are each a product of 
operators of the form Wiiji^r, or Wi'j'^r, where r stands for a time variable in {ti, . . . ,tk,ti, . . . ,tk'}. 
Moreover, the product ABC contains k W^s and k' Ws. Finally, each time variable in 
ti, . . . ,tk,ti, . . . ,tk' appears exactly once in the product ABC. 
Let if e n'^iP+''+''') and estimate 



E E *s*2' 

CG.S(p,fc,0) S'eJ2(p+/c,fc',0) 

^ E , , , 



dt 

A'x(t) Jo 



dt / ds 



ds J^^^ dt! Aa[P!p^, BP!F+'>^,{s)W,j,,Cip 
dt' 



[0,*]*=' 



Aa^ll^,BP[^+''^,{s)W.,,,C^ 



We now perform all time integrations, starting from the left, and using at each step the estimate 
as well as 

ft 



dr \\Wi'j' ^r^^W ^ \/2'KK'^t ||(p| 







which follows trivially from ([5^ . Also, Lemma f2. II implies that |l7(s)|l = ||7|1. Thus we find 
that 



^< E E 

Qei2(p,/c,0) Q'G^(p+fc,fc',0) 



(2^.^2t)''/'||a('')||||7||||(/^l! 
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Using the bound 



which can be inferred from ()29p . we find 

^2^\ {k+l)/2 



1&^/2ttkH {QA^2TinHf {?,2V2^tf' Haf")]! ||7|| 



Let t < (2"7rK2)-i. Now Lemma O implies that ||^®(p+/c+A;')j](p+fe) ||^ <^ i Ugj^g the inequal- 
ity Tr(^r) ^ ||^||||r||i we therefore find that (gl]) is bounded by 

k' , . (32^2^)'' , , 

ii„(p)iiiu,ii _ i«P V V ; iia(f)|l|l7|l . 



16P ^(32^2^)'' (l6\/2^)'' lla^P^llhll = 16^-!: 



l&^/2'KnH 



The second term of Rij{t) is equal to the complex conjugate of the first. We thus arrive at 
the desired bound 

\RUt)\ ^ 2.16^i^^^^^^ll||a(^)||||7||. (43) 
Therefore the last line of ((in|) is bounded by 

oo 



1 - 1Q^/2ttkH ^.^p 



1 - 1&^/2ttkH (1 - 32V2^^; 



where we used the estimate Yll^=o{P + k)^ x'' ^ (Xz^yr+T- 

Next, we note that the second line of (HO]) , i.e. the rest term, vanishes in the limit K ^ oo. 
The procedure is almost identical to (in fact easier than) the above estimation of |i?*j (t)|. The 
result is 



dt Tr (C^f "Ha(f))7(tK)®'''+'^^sL^+^) 

as K ^ OO. 

Summarising, we have proven: 

Lemma 5.2. Let a^^^ e B(t6]^^). Then, for small times, 



{i2V2^t) „ , 
1 - IQ^/2tthH 



k=0 



^ \\a^P^\\MC{p,K,t), 



for some constant C{p, k, t). 



5.3 The mean-field limit 

We now have all the necessary ingredients to prove our main result. Take some fixed orthonor- 
mal sequence $ = (iy9i)igN- Denote by the truncated sequences = i^i, ■ ■ • , Vn), and let 
$Ar(i) be the solution of the Hartree-Fock equation (|20p with initial data $Ar. The iV-particle 
density matrix evolved with the Hartree-Fock dynamics is 

fjv(t) := |5($Ar(t)))(5($^(t))|. 

The A'^-particle density matrix evolved with the microscopic dynamics is 

T^it) e-'*«~|5($jv))(5(<i>Ar)|e"^". 
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The p-particle marginals are defined by 

rj^^(t) := Trp+i...^rjvW, 
The one-particle density matrix satisfying (113p is 



The quantities T^'{t) = I]^^7Ar(t)®P and r^^(i) are asymptotically equal: (EH) implies that 

„2 



oil. * ^- 



(44) 



Next, let a(P) G B{r6^^). From dMl) and ([331) we get, for small times t 
Tr(Ajv(a(^'))rAr(i) 



OC 

oo 



Using Lemma [521 and (EH) we therefore get that, for small times, 



Tr 



(ai^)7^^)(0)-Tr(A^(a(-))r^(t))| < ^^P^\\a^^^\. 



Since the quantum-mechanical and the Hartree-Fock time-evolutions preserve the trace norm, 
we may iterate the above result, like in [9], to get: For all times t G M we have that 

Tr(a(P)7j*')(i)) -Tr(A^(aW)r^(t)) I < C{p,K,t)\\a^^^\\ f{N) , 

with limjv^oo / (A^) = 0. Thus, by using the duality B = [L^)* we find 



Theorem 5.3. Let p e N and t e M. Then 

\\V^^^it)-T^S\t)\l 
In particular, if a^^^ G B{l-6^^), we have 



0, 



as N ^ oo. 

Remarks. 1. The limit A^ ^ oo of r^''(t) does not exist in ||-||i. Indeed, limAr_+oo||r^''(i)|| = 
but Trrj^^(t) = 1 (similarly for f^^(t)). 

2. As in [9], one can show that the function / is a power law: f{N) ~ N~^^*\ with (3{t) > 
for all t. However, P{t) ^ as t ^ oo. Our bound on the rate of convergence is therefore far 
from the expected optimal rate P{t) = 1, which we only obtain for short times. 

3. Although the exchange term — * {fi'Pj))'Pj is essential for our proof, it is not 
clear from our analysis whether it is of leading order as A^ — > oo. The exchange term is known 
to be of subleading order in the scaling of [14], and hence in that case it does not play a role in 
the limiting dynamics (see [6]). 



6 A Egorov-type result for small times 

In this section we describe how the many-body dynamics of fermions may be seen as the quan- 
tisation of a classical "superhamiltonian" system, whose dynamics is approximately described 
by the Hartree-Fock equation. 
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6.1 A graded algebra of observables 



We start by defining a Grassmann algebra of anticommuting variables over the one-particle 
space Ti. — L^(M'^), and equip it with a suitable norm. Formally, we consider the infinite- 
dimensional Grassmann algebra generated by {^{x), ipi^)}xeR'' ■ it turns out, this algebra can 
be made into a Banach algebra under a natural choice of norm. This norm is most conveniently 
formulated by identifying elements of the Grassmann algebra with bounded operators between 
L^-spaces. 
Let 

a = (a^^^'^WN, a(P'«) eS(HL'\7^L^^ (45) 

be a family of bounded operators. Such objects will play the role of observables in the following. 
By a slight abuse of notation we identify a'-^^'?-' with the family obtained by adjoining zeroes to 
it. 

Define 

'B'^ {a = (a^P^^)) . ^(p^'?) = for all but finitely many (p, q)} . 
We introduce a norm on Q5'~^ through 

\\a\\^a := ^|la(^-'^)||, (46) 

and define S'' as the completion of 03*^. 

We also introduce a multiplication on ^B*^ defined by 

(a6)(P'«' (-1)^^(^1+91) P_(a(Pi'«i) «)6(P=^'«=^))P_ . (47) 

Pl+P2=P 

91+92=9 

The seemingly odd choice of sign will soon become clear. It is now easy to check that QS*^ is an 
associative Banach algebra with identity 

Note that QS*^ bears a Z-grading, with degree map 

dega^P'''-' := p — q . 

An observable is gauge invariant when its degree is equal to 0. One readily sees that 

ab^ (-l)'i°s«d<=g65a. 

We now identify with a Grassmann algebra of anticommuting variables. For f ^ H 
define V(/) e B{H; C) C W through 

W)9 ■■= if, 9) (48) 

and i^if) e B{C- H) C »^ through 

^(/)z := fz. (49) 

We may now consider arbitrary polynomials in the variables {ip{f),ip{f) ■ f € Ti,}. It is a 
simple matter to check that 

V'(/)V^(ff) + ^(.9)V'(/) = V'(m.9)+?(5)V^(/) = Hf)JK9) +i^{9)JKf) = 0, 
for all /, g G Ti.. Furthermore, we have that 

i^ifp) ■ --Hfimgi) ■ ■ ■ ^(59) = P-'^ 1/1 ® • • • ® ® • • • ® 59I ^-'^ ■ (50) 

Linear combinations of expressions of the form ([50]) are dense in S*^. It is often convenient to 
write a family a of bounded operators using the "Grassmann generators" {?/;,■(/'}• To this end 
we set 

ip{x) := -0(4), ^(x) := tp{Sx), 
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where Sx is Dirac's delta mass at x. Expressions of the form (j50|) are now understood as densely 
defined quadratic forms. One immediately finds 



- aG 



X i^{xp) ■ ■ ■ iIj{xi) a^P'''\xi, ...,Xp;yi,..., yq) tl){yi) ■ ■ ■ il^iyg) . (51) 

We use the notation AP{a) to emphasize that the family a is represented using Grassmann 
generators. 

6.2 A graded Poisson bracket 

Next, we note that *B'-^ carries the graded Poisson bracket 



{a,b} 



i / da; 



6 



S , 6 

■b + a- 



6 



Si}{x) 6ijj{x) 6i;{x) Si}{x) 



(52) 



where a,b £ S*^. Here we use the usual conventions for derivatives with respect to Grassmann 
variables (see e.g. [17], Appendix B). In terms of kernels the graded Poisson bracket can be 
expressed as 



{V'(a;),V'(y)} = iS{x-y) {ip{x),t(;{y)} = {V'(a;), ?A(2/)} = 0. 
We now list the important properties of the graded Poisson bracket. 

(i) {a,b} = (-l)i+d'=s'^ 'i°s''{&,a}. 

(ii) (-l)'^'=Sb(dcga+dcgc)|^ |^^^||_|_(,y(,li^ pgj.j^^^j^^jQjjg ^ 

(iii) {a, be} = {a,b}c + {-l)'^''ii'"^''^''b{a, c} . 
Proof. Let us start with (i): 

Sa Sb 



(53) 



{a, b} — i J dx 



i / da; 



f j^Ndcgg+l , / j^Ndoga+l 

^ ' diPix)S^{x) ^ ^ Si^{x)6t^{x) 

^_2^dcga dcg f)+dcg b ^ ^_-[^-jdoga dcgb+dcgf) ^ 



5^{x) 5tP{x) 



6i){x) 5il:{x) 



= (-1) 



1+dcg a dcg b 



{b,a}. 



In order to show (ii), we note that the left-hand side can be written as a sum of three terms, 
the first of which contains second derivatives of a, the second second derivatives of b and the 
third second derivatives of c. Let us consider the third one. It is equal to the terms containing 
second derivatives of c of 

(_;L)dcg6(doga+dcgc)|^^ _^ ^_-^)dcgc(dcgb+dcga)|^^ ^|| 

^ (_;^)dog6(doga+dcgc)|^^ ^ (■_-^^dcgc(dogb+dcga) + l+dogc doga|^^ ^ 

where (i) was used. Define the derivation Lab := {a, 6} . Thus we need to compute the terms 
containing second derivatives of c of 



dcg a dcg6+dcg6 dcgc 



LaLbC- (-1) 



dcg b dcg c 



LbLaC . 



Since we are only considering terms containing second derivatives of c, both derivations La and 
Lb must act only on c, and one finds 



( — 1)'^°^'^ dcg6+dcg& dcgc^^^^^^ _ ^_-|^'jdcga dcg 6+dcg 6 dogc^^^^^^ 

We omit the straightforward proof of (iii). 



0. 



□ 
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Furthermore, one finds by explicit calculation 

= i (-l)(P^+i)(Pi+«i)gip2 [(a(Pi''?i) (g) (1(91-1) ^ 5(P2,«2)) 

- i (-l)(«i+i)(''^+''^)pig2 ^ i(pi-i)) (1(92-1) ^ a(Pi,9i)^l , (54) 

6.3 States 

To the algebra of observables (25"-^, IHIsbg) is associated the space of states 

$n := (s«)*. 

By using the standard argument of the proof that (/^)* — l°° (see e.g. [15]), one finds that 
9^ = {p = {p^P^'^X^.^n : gS(wL'^hL^^), ||p|k<oo}, 

where 

IIpII^h := sup||p(f'«)||i 

p,q£N 

and 

sup{|Tr(/f'«)a(«'P))| : a^'-f^ G hL'^) , ||a(«'P) || ^ l} . 

Note that if p = q then ||-||i is the usual trace norm. The dual action is given by 

(p,a) J2 Tr(p(P''?)a(«'P)). 

p,gGN 

We abbreviate p(P'P) = ^(p) in the case of gauge invariant states. Next, we note that ([50]) implies 
that the operator kernel of p(P'9) may be expressed as 

p^P'''\xi,...,Xp;yi,...,yq) = (p ,il!{yg) ■ ■ ■ ■il;{yi)ip{xi) ■ ■ ■ il;{xp)) . (55) 

There is a particular subset of gauge invariant states that is of interest for studying the 
Hartree-Fock dynamics. Let 7 € I? be a one-particle density matrix. Define the state p^ 
through p^'"^^ — if p ^ q and 

p\P^ := 7^^^ , (56) 
where 7^^^ is defined in (fTT]) . One immediately finds Hp-ylji — 

6.4 Hamilton function and dynamics 

Let h be the one-particle Hamiltonian and w the two-body interaction potential. We define a 
Hamilton function on (a dense subset of) the phase space D\ through 

H := A^{h) + ^A^{w) 

dxdy ^p{x)h{x;y)^/j{y) + ^ [ dx dy i{;{y)ip{x) w{w - y) ^p{x)^p{y) . (57) 



2 _ 

The Hamiltonian equations of motion read 

a = {H,a}, 

where a S S*^. Instead of the "Heisenberg" evolution of a we consider the dual "Schrodinger" 
evolution of states: 

(pit), a) {p,ait)). 
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The equation of motion for states reads 

idtp'-P'''\xi,...,Xp;yi,...,yg) = i^h^i -^hy^ 1 P^^'"^^ (a^i, • • • , a;p; yi, . . . , y^) 

+ du i ^ w(u - Xi) - ^ w(u - y^) . . .,Xp,u;yi, 

•' \i=l i=i / 



yq,u). (58) 



This has the form of an infinite hierarchy, which decouples over subspaces of different degree. 
In order to show ()58p we compute 



i{7J, i:(yq) ■ ■ ■ ip{yi)ip{xi) ■ ■ ■ ip{xp)} = i^^K^ -"^hyA ip{yq) ■ ■ ■ ij;(y i)ip{xi) ■ ■ ■ ip{xp) 

+ / ^" I X! "'(^ ~ ~ X! ~ I V'(w)V'(y</) • • ■ i^{yi)i^{xi) ■ ■ ■ il){xp)i}){u) . 

\i=l i=l I 



Then ([581) follows from ([55]) and 

iatp(^'''?^(xi, . . . ,Xp;yi, . . . = idt{p,ip{yq) ■ ■ ■ ip{yi)il;{xi) ■ ■■ip{xp)) 

= {pA{H,tp{yq) ■ ■ ■tp{yi)ij{xi) ■ • •?/'(xp)}) . 
Next, we outline how to solve the equation of motion (|58p . Let us first rewrite it as 

p 9 



i=l 1=1 

+ f]Trp+i^,+i(iy,p+ip(^+i.^+i)) - ^Trp+i,,+i(p(P+i^«+i)W^,,+i) . 



We may now proceed exactly as with the density matrix Hartree-Fock equation p2p . i.e. express 
it as an integral equation in the interaction picture. This yields a tree expansion for the quantity 
Xr(p(P'9)(i) a^'^'P)), where p(0) e 91. We omit the uninteresting details. As above, the tree 
expansion converges if t < T, where 

T := (2"^k2)-i_ (59) 

Unfortunately, the time evolution (|58p does not preserve the norm of p, which means that we 
cannot iterate the short-time result. 

From now on, we only consider gauge invariant invariant quantities. Take some gauge 
invariant state p — {p^P'>)p^fi S 91. For simplicity, we assume that the sequence p is finite (as is 
the case if p is defined by a Slater determinant, see below). Let us denote the Hamiltonian fiow 
on by We have seen that (/)* is well-defined by its tree expansion for t < T. The solution 
of ([55]) with initial data p, p{t) = (?!)*(/9), satisfies the equation 

p(p){t) = p^P>-i I ds > 'Trp+i|iy,p+i,„p^f+^J(.s)| , (60) 



(t) = P^P^-i I ds ^Trp+i[iy,p+i,„p(f+i)(.s)] , 
■^0 i=i 



where p^P^{t) :— e'^> p(p)(t) e '^Zi'^i*, Let us take a gauge invariant observable g'^P'P^ 
a(p) e Sf^, where 

{a e : a^P'^^ ^ if p ^ q} 
is the set of gauge invariant observables. Then (j60|) implies 

TT{a<-P^p^P^t)) = Tr(aiP^p(f)(t)) 

= Tr(4Vp))-|-i fds VTr([M/.p+i.„4^)®llp(P+i'(.)). 

•^0 *=i ^ ^ 
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Iteration of this identity gives 

oo 
k=0 

Summarising: 

= f;Tr(G^"^(a(P))/P+'=)) 

k=0 

This series converges for t <T, uniformly for bounded ||a^^''||!gG and Therefore we get 

the norm-convergent series 

AG(a(rf)o0* = f;A«(G^")(a(j'))), (61) 

k=0 

provided that t <T. 

Finally, we discuss the relationship between the Hartree-Fock dynamics and the dynamics 
generated by (j58p . Take a density matrix 7 S I? and consider the state p — p^y defined in 
([SS]). If one chooses a sequence such that ||7Ar|| ^ as TV — > 00 (e.g. a sequence of Slater 
determinants), then Lemma [5.21 implies that ([55|) and the Hartree-Fock equation describe the 
same dynamics for large N. 



6.5 Quantisation and a Egorov-type theorem 

In this final section we introduce a Wick quantisation of the above "superhamiltonian" system 
and formulate the mean-field limit as a Egorov-type theorem. From now on, we use n to denote 
the number of particles, and N to denote a free parameter (the inverse "deformation parameter" 
of the quantisation). Ultimately, we shall choose n = N. The underlying ideas were described 
in detail in [9]; here we merely state how they apply in the current setting. 
Consider the fermionic Fock space 



(») 



on which act the creation and annihilation operators a*{x),a{x) (see e.g. [4]). We define the 
rescaled operators as ip'^ix) := a'^{x)/\/N, where # € {*, nothing}. They satisfy the canonical 



anticommutation relations 



ip*^{x),i!%{y) 



0, 



1 



j:^^{x-y), 



where [A, 5]+ := AB -f BA is the anticommutator. 

Let a(P) e B{T-6p^) with distribution kernel a^P\xi, 
sation of a(P) is defined by 



^iv(a;),i/'Ar(2/)J^ 
. . , Xp] yi, . . . , yp). The second quanti- 



A^(a(P') := J dxi---dxpdyi---dyp 

'$*n{xp) ■ ■ ■ ^*n{xi) a^''\xi, ...,Xp;yi,..., yp) ?/'7v(2/i) • ■ • ^/'7v(2/p) ■ 

This defines a closed operator on Moreover, restricted to the subspace Ti.'"_^\ A'^{a'^P'>) is 
equal to Ajv(a(^')) defined in (^5]) . Additional properties of the operation A^ are listed in [9]. 
Let Sl*^ be the linear hull of 



{A^(a(f)) : a(P) eB'^P\peN} 
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We define quantisation as the linear map : St*^ — > 21*^ defined by the formal replacement 
ipix) '—^ 'iPn{x) and ip{x) i— *■ il)*j^{x) followed by Wick ordering. In other words, 

V)^ : A«(aW) - A«(a(^')). 



Using it is easy to see that, for A,B^ 21*^ 

- - AT-i 



{A,B}^ + 0{N-^). 



This identifies TV -'^ as the deformation parameter of (Oat- 

Extending the definition of (•)^ to unbounded operators in the obvious way, we define a 
Hamiltonian Hjq on as the quantisation of the Hamilton function H defined in (|57p . When 
restricted to 7i^\ NHn is equal to the Hamiltonian with mean- field scaling 

Now , (|34p and (pTjl yield the following Egorov-type theorem. 

Theorem 6.1. Let A e 2t<^ and t < T, with T defined in Then 



\tNHN 



Ang 



-UNHn 



{Ao(j)t 



N 



C_ 



for some C > 0. 



A Hamiltonian formulation for density matrices 

In Section [21 we chose a Hamiltonian formulation of the Hartree-Fock equation ([7]) in terms of 
sequences of orbitals. Alternatively, we could just as well have used a Hamiltonian formulation 
in terms of density matrices. To see how the density matrix Hartree-Fock equation (fT2|) can be 
written as a Hamiltonian equation of motion of a classical Hamiltonian system, consider the 
Hilbert space 

the space of Hilbert-Schmidt operators, with scalar product 

:= Tt{k*p). 

We write the density matrix 7 e I? as 7 = kk*, where k G 7i. The classical phase space is then 
given by a Sobolev-type space of Hilbert-Schmidt operators 

r := [neU : Tr(K*(l - A)k) < 00] . 
We define polynomial functions on F through 

B(a(P))(K) := (K®P,a(P)K®P), 

where a^P^ G B{n®P). 

The affine space F carries a Symplectic form defined by 

UJ = —i J dx dy dR,{x , y) A dK{x , y) , 

where K{x,y) E L^(M'^ x M.^) is the operator kernel of k. The Poisson bracket then reads 

{K.*{x,y),K*{x',y')} = 

{K{x,y),K{x',y')} = ~iS{x ~ x')S{y -y') . 

The Hamilton function is defined by 

H := B(/i) + ^B{W) . 
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By using Sobolev-type inequalities one readily sees that H is well-defined on F. After a short 
computation, one finds that the Hamiltonian equation of motion, 

SH 

idtK{x,y) = — = i{H,K{x,y)}, 

dK{x,y) 

reads 

idtK = /iK + Tr2(yVK(g) (kk*)) . 
It follows that 7 = KK* satisfies 
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